Complex oscillatory yielding of model hard sphere glasses 
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Abstract 

The yielding behaviour of hard sphere glasses under large amplitude oscillatory shear has been 
studied by probing the interplay of Brownian motion and shear-induced diffusion at varying os- 
cillation frequencies. Stress, structure and dynamics are followed by experimental rheology and 
Browian Dynamics simulations. Brownian motion assisted cage escape dominates at low frequencies 
while escape through shear-induced collisions at high ones, both related with a yielding peak in G" . 
At intermediate frequencies a novel, for HS glasses, double peak in G" is revealed reflecting both 
mechanisms. At high frequencies and strain amplitudes a persistent structural anisotropy causes 
a stress drop within the cycle after strain reversal, while higher stress harmonics are minimized at 
certain strain amplitudes indicating an apparent harmonic response. 

PACS numbers: 
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Hard Sphere (HS) colloids have been used as model systems to study a plethora of funda- 
mental condensed matter problems such as the interplay between equilibrium phases (crys- 
tals and liquids) and non-ergodic states (glasses or gels), and their behavior under external 
fields such as shear A major goal is to develop an understanding that will enable 

tailoring of the mechanical and flow properties based on the structure and dynamics at the 
particle level. HSs form metastable glasses above a volume fraction of about ~ 0.59 where 
crystallization and long-time diffusion are suppressedH , y] . Such states exhibit solid-like be- 



havior at low stresses/strains and shear-melting (yielding) above a yield stress/strain [5|. 
This phenomenon is investigated more effectively with a combination of rheological and op- 
tical or scattering techniques that may unravel the link between microstructure, dynamics 

, IB| • Such studies have shown 

that particle cages are deformed under shear and particles move irreversibly when a critical 
strain is exceeded, leading to flow. 

Experimentally, large amplitude oscillatory shear tests are widely used to monitor yield- 
ing due to their simplicity and relation with linear elastic and viscous moduli, G' and G" . 
However, their interpretation in the non-linear regime becomes complex due to stress distor- 



tion introducing higher harmonics 



10|. In a wide range of systems (emulsions, polymers and 



colloids) a generic peak in G'l (viscous modulus of the fundamental frequency) is observed 
representing an increased ener gy d issipation near the yield strain, 7^, where G'^ = G'(, be- 
yond which the sample flows jolJllll. On the theory side. Mode Coup ling Theory (MCT) jl2| 
as well as the semi-phenomenological Soft Glassy Rheology (SGR) [13| are able to capture 
some aspects of Large Amplitude Oscillatory Shear (LAOS) such as the peak of G". Nev- 
ertheless, the underlying mechanisms relating stress with shear induced structure such as 
cage deformation, breaking and reformation, as well as particle displacements, are poorly 
understood, and the frequency dependence is largely unexplored. 

LAOS experiments are expected to provide valuable information on the interplay of Brow- 
nian motion and shear during yielding of HS glasses. Along this line light scattering-echo 
experiments probing the average particle displacements under oscillatory shear 8|, |9j revealed 
a transition to irreversibility beyond a critical strain amplitude js, This can be viewed 
as the analogue of the shear induced irreversibility observed in concentrated non-Brownian 
particles {14 1 . since in colloidal glasses out-of-cage diffusion is frozen, while in non-Brownian 
particles diffusion is absent at all length scales. While experiments combining rheometry 
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with scattering or microscopy are rather demanding [8|, Il5l-ll8|. computer simulations can 
provide an alternative route. 

Here we use a combination of oscillatory shear rheometry and Brownian Dynamics 
(BD) simulations to investigate the links between structure and particle dynamics with 
the non-linear rheological response of HS glasses in a wide range of frequencies, u, non- 
dimensionalized by Pe° = utb, with tb = R^/Dq, R the radius and Dothe free diffusion 
coefficient, or Pe^ if the short-time self diffusion coefficient Ds{ip) is used. At low Pe° yield- 
ing is related to Brownian-assisted irreversible particle motion manifested in a Dynamic 
Strain Sweep (DSS) with the peak of G'/. However, at the largely unexplored regime of high 
PeJi, we detect collision-dominated yielding and, within the oscillation period, a strongly 
anisotropic structure causing a reduced stress beyond strain reversal (memory of structure 
due to lack of Brownian relaxation). At intermediate Pe°, the sample is affected by both 
mechanisms as manifested by a novel, for HS glasses, double peak in G'/. 

We used sterically stabilized poly(methyl methacrylate) (PMMA) model nearly hard- 
sphere particles with two radii, R = 358 nm and 130 nm (polydispersity 10 — 12%) suspended 
in octadecene and octadecene/bromonapthalene mixture respectively, in order to expand the 
Pe^j range. We prepared different volume fractions, by diluting a random close packed 
)atch with the exact (f then adjusted by matching G[ of the two systems in agreement with 
lol . Experiments were performed on an ARES strain controlled rheometer (with 25mm 
diameter /O.Olr ad angle cone-plate) and a solvent trap to eliminate evaporation. In BD sim- 
ulations, HS interactions were implemented through the potential free algorithm 2o| . Oscil- 
latory shear was applied with periodic boundary conditions using typically 5405 particles 
with 10% polydispersity to avoid crystallization. 

In figure 1 we show DSS tests performed at low and high Pe^^, together with the evolution 
of the normalized intensity of all higher stress harmonics, laii/h = '^li/Ii, = 2n + 1, 
n > 1) , as well as Lissajous curves (intracycle stress versus strain) for representative strain 



amplitudes, 79. The low Pe^^ (= 0.5) data (fig. la) exhibit the typica. 



DSS response with 



a G'{ peak around 7^^ (at G'^ = G'{) similar to previous studies |lll. llOL |21| . The non-linear 
response is accompanied by progressively larger intracycle non-linearities as indicated both 
by the Lissajous curves and the increasing laii/h as expected |10l. Il9| . The latter increases 
beyond yielding and reaches almost 30% at high 7q as found previously jl9|. Moreover, the 
Lissajous plots show a transition from a linear viscoelastic behavior (elliptical shape) at low 
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FIG. 1: Top: Dynamic strain sweeps for (a) R = 130 nm, ip = 0.60 at w = Irad/s (-Pe^, = 0.04, 
Pe^j = 0.5 ) and (b) R = 358 nm ip = 0.60 at a; = 1 rad/s (Pe° = 0.9, Pe^; = 11.2), with the 
1st harmonic of the elastic G'l and viscous G" modulus as a function of strain amplitude. Middle: 
Representative Lissajous plots are shown in different strain amplitudes as indicated. Bottom: 
Normalized total intensity of the higher harmonics of the stress, laii/Ii- 



7o, to a paralle 



at 7o > 7, 
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ogram pattern indicative of a intracycle sequence of elastic-plastic response 



19|. 



At high Pe^)-, achieved with large particles {R = 358 nm) at a; = 1 rad/s [Pe^^ = 11.2) the 
response is qualitatively different (fig. lb). Firstly, the peak of G'{ shifts to higher 7q, beyond 
7y. Secondly, laii/h (and individual hn+i/h) exhibits a non-monotonic behavior showing 
a first maximum around ■jy and subsequently decreases substantially well inside the non- 
linear regime. Hence, the sample exhibits a more harmonic stress response (anharmonicity 
is lowered) even though under non-linear LAOS. 

We further explore the Pe^^ dependence by changing u while keeping 7q constant. In 
fig. 2a experiments with varying u (at 7q = 100%) reveal the transition from the low Pe^^ 
rectangular shaped Lissajous curves refiecting a sequence of elastic and plastic responses, to 
the high Pe^j regime with a characteristic ellipsoid with a double concave distortion caused 
by reduced stress in the quandrants II and IV after strain reversal. The intensity of the 3rd 



harmonic at 100% (above 7^^) exhibits a minimum with u in experiments, = 0.62, and 
BD, if = 0.60, (fig. 2c); note that the position of minimum is if and 7q dependent. While 
in both Pcui regimes the stress response is highly anharmonic, with significant /3//1, during 
the transition the Lissajous curves (fig. 2a) acquire an ellipsoid shape involving almost 
zero higher harmonic contributions. BD simulations showing identical rheological response 
(fig. 2b) with experiments are able to provide valuable structural information revealing 
the underlying mechanism of such stress reduction. In fig. 2d we plot the 2D projection 
of the pair correlation function in the velocity-gradient (xy) direction, g^ylr), at specific 
points inside the oscillation cycle for Pe^j = 100, similar to findings under steady shear j?]. 
Contrary to what we find at low Pe^j (supplemental material), here the structure is highly 
anisotropic at the point of maximum strain (zero shear rate; point A in fig. 2b). Moreover, 
such anisotropy is persistent during a large part of the successive quandrant where the shear 
has been reversed (points B and C). Such an anisotropic cage, created during high Pe^ 
shear in one direction, allows fiow with less stress (due to fewer particle collisions) when 
shear is reversed. Thus, in quandrants II and IV the stress is reduced if compared to a fully 
harmonic viscous response corresponding to the fiowing anisotropic structure of quadrants 
I and III (dash-dot line in fig. 2b). Only beyond zero strain (points D and E) is the 
structure reversed and the stress comes back to the maximum values within the period (see 
supplemental material). Such response is absent at low Pe^ since Brownian motion relaxes 
shear-induced structural anisotropy more efficiently, and the stress response is similar in all 
quandrants. 

The data presented above verify the existence of the two Pe^ regimes: The low Pe^^ 
one, conventionally studied up to now, where Brownian motion is dominant and the high 
Pe^j where shear-induced particle collisions introduce novel LAOS features related to the 
persistent structural anisotropy and the consequent reduced stress after strain reversal. The 
transition from Brownian activated yielding, where particles under shear escape their cages 
assisted by thermal motion, to collision-induced cage breaking at high frequencies is linked 
to pronounced irreversible particle rearrangements and decreasing 7^ at low frequencies [sl,^]. 

We further investigated HS glasses at higher if and frequencies corresponding to an inter- 
mediate {(f dependent) Pe^^ regime. Experimental DSS, Lissajous plots and laii/h shown 
in fig. 3 as a function of 7o for if = 0.639 reveal an even richer mechanical response at 
Pe° = 0.4 [Pci^ = 8). A main observation here is the unambiguous detection of two peaks 



5 



FIG. 2: Lissajous plots at 100% strain amplitude (a) for experiments at different Pci^ denoted by 
the dash lines in (c) with the 4 quantrants of the oscillatory cycle indicated and (b) BD simulations 
at Pe^j = 100. The dashed line in (b) indicates a viscous harmonic stress strain response (see also 
supplemental material), (c) Pci^ dependence of the 3rd harmonic at 100%, for experiments at 
ip = 0.62 (R = 358 nm) and BD simulations at ip = 0.60. (d) 2D projections in the velocity- 
gradient (xy) plane of the difference of gxyif) under shear from that at rest from BD simulations 
at indicated points within the cycle, points A-E in (b). 



in G'l detected for the first time in HS glasses. Such a feature was so far observed only in 
attractive glasses and gels indicating a two-step yielding due to two length-scales present in 
attractive systems i.e. the interparticle bond and the cage or cluster size 2l|, |22|. However, 
the double G'( peak seen here must be of a different nature since a second length-scale is 
absent and moreover, the phenomenon is only observed in a narrow range of Pe^. The first 
peak of G'l is identified with the one observed at low Pe^^ by a direct comparison of the two 
DSSs (-Pe^ = 0.04 and 0.4) shown in fig 3a. The two peaks signify the maximum in energy 



FIG. 3: (a) Dynamic strain sweeps with R = 130 nm particles at = 0.639 and an intermediate 
= 10 rad/s, Pe° = 0.4, Pe^ = 8 thick black lines) and low = 1 rad/s, Pe° = 0.04, Pe^ = 0.8 
thin red lines) frequency regime, (b) Indicative Lissajous plots for Pe^^ = 8 at different strain 
amplitudes and (c) Normalized total intensity of higher harmonics laii / h ,versus strain amplitude 
for Pe^ = 8. 

dissipation during the two yielding mechanisms at low and high Pe°, attributed to cage 
breaking via shear- assisted activated hopping and through particle collisions, respectively. 
The Lissajous figures and higher harmonics reveal a transition from the low 7q linear re- 
sponse to a viscoplastic fiow at high 'Jq passing through two states with apparent harmonic 
response as indicated by the two minima of laii/h- At these strain amplitudes the Lissajous 
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FIG. 4: BD simulations at ip = 0.60: (a) DSS data showing G'l and G" as a function of 7q for 
different Pe^j as indicated, (b) Power law exponent for the Jq dependence of and D(,ff{t = T) 
as a function of Pe^) (c) Average D(,ff(t = T) versus 7q for different Pe^^ and (d) D(.ff{t = T) 
versus Pe^) for 7q = 10%, 30% and 100%. The corresponding Deff{t = T) for a system at rest is 
also indicated. 



curves acquire an nearly ellipsoidal shape due to the compensation of the structural phe- 
nomena during the transition from low to high Pe^) (as in fig 2c). Note that for this sample 
{R = 130 nm, (f = 0.639) the high Pe^^ (fig lb) was not within the experimental window. 

Fig. 4a shows G[ and G'l from BD LAOS tests for ip = 0.60 at different Pe^. The 
dependence of G[ and G" at high Jn{> ly) follows a power law decrease G[{G'()o: 7^' 
(7g") as detected experimentally 111. Il9|. While Maxwell-type models give u' = Iv" = —2 
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and MOT (around the glass transition) predicts lower values but similar i/' /y" ratio 
experiments in HS glasses show deviations from such simple dependency 19|. In agreement 
with experiments (fig. lb and supplemental material), BD simulations give Pe^j depen- 
dent exponents (fig 4a) with u" approaching —1 and at low and high Pe^^ respectively. 
Therefore, at high Pe^, where collision activated out-of-cage particle rearrangements are 
dominant, 7oG'/ ^ cr(7max)) ^ measure of energy dissipation per unit strain, is proportional 
to shear rate {'JoUj) similar to the limiting high shear rate viscosity behavior under steady 
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shear. On the other hand at low Pe^j a G'( oc 7q ^ dependence corresponds to a steady shear 

. -1 

shear thinning response of a HS glass at the yield stress plateau with 77^^^ oc 7 . Note 
that at this stage, the absence of a double G'l peak in BD can not be firmly attributed to 
the absence of hydrodynamic interactions, since BD at high volume fractions could not be 
conducted at a similar number of data points as experiments due to computational time 
restrictions. 

Further insight in the two yielding mechanisms is gained by examining microscopic par- 
ticle dynamics within a LAOS cycle by BD. In Fig. 4c we show the effective diffusivity, 
Deffit = T) = (A2^(T))/T with {Az'^(T)) the mean square displacement in the vorticity 
direction, z, (at t = T ), as a function of 7o for several Pe^. At low Pe^, Df.ff{T) increases 
sublinear with 7q, whereas as Pe^ is increased it exhibits progressively a weaker increase at 
small 7o and a stronger one at higher. The constant Deff{T) at low 7g indicates that prior 
to yielding in-cage diffusion is unaffected by shear. For 7q > 7^, Deff(T), corresponding 
to out-of cage diffusion, increases sublinear (7q = 100%, fig. 4d) with a power law expo- 
nent ^ 0.8 at low Pe^j, which approaches 1 for Pe^ > 1 where direct particle collisions 
are dominant. In comparison, the corresponding Deff{T) for 7q = 10% and 30% show 
an initial sublinear dependence at low Pe^ and subsequently approach the curve at rest, 
since for such frequencies 7^ > 30% [9I and the sample has not yet yielded. Note that a 
similar power law increase (with exponent ~ 0.8) has been detected in HS glasses at low 



steady shear rates by confocal microscopy 



a], a weak but systematic deviation from the 



linear MOT prediction 2J| and closer to agreement with non-linear Langevin equation the- 



ory involving activated hoping mechanisms 



25| . For our oscillatory BD a linear dependence 



26|. 



is reached for Pci^ > 10 as expected for non-Brownian particles under steady shear 
We suggest that the sublinear and linear dependencies probed here at low and high Pe^j, 
respectively, reflect the two different mechanisms involved in the two regimes. Note that 
similarly with -De//, 7o^i increases linearly with 7q at high Pe^; the linear 7q dependence 
of both quantities is linked to collision-induced yielding. On the other hand at low Pe^ both 
quantities increase sublinearly, -De// oc 7g* and JoG'l oc ■y^-'^ (the latter exponent tends to 
zero as Pe^ is lowered) reflecting plastic flow, with a power law stress behavior, and complex 
Brownian/shear-activated particle hoping. 

In summary, the combination of experimental oscillatory rheology and BD simulations 
has revealed the complete mechanical fingerprint of HS glasses and the related underlying 
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microscopic structure and dynamics over a wide frequency regime where both Brownian and 
non-Brownian behavior is probed. At low Pe,^, commonly studied up to now, Brownian- 
assisted irreversible motion takes place during yielding with a single peak of G'l and strong 
higher harmonics of the stress at large 7q. In this regime the microstructurc under shear 
is only weakly anisotropic while the shear-induced diffusivity scales sublinear with Pe,^. At 
high Pe^^ yielding is dictated by collision-induced displacements linked with a shear-induced 
long time diffusion that increases hnearly with Pe^. A single G'[ peak is detected at 7q > 
which eventually turns into a plateau at the limit of large Pe^^. The structure under shear is 
strongly anisotropic but more interestingly exhibits a hysteresis under strain reversal. Such 
structural memory due to lack of Brownian relaxation, causes a characteristic stress drop 
after strain reversal as revealed by the Lissajous curves. At some characteristic 79, higher 
harmonics drop to almost zero indicating an unexpected harmonic response even though the 
sample is under non-linear shear. Finally, at a (/? dependent intermediate Pe^^ regime the 
sample is affected by both mechanisms (Brownian- and collision-induced yielding) resulting 
in a double peak of G'{ and two minima of higher harmonics, a feature that is detected for 
the first time in simple hard sphere glasses. 

The rich mechanical response of a model hard sphere glass revealed here in oscillatory 
shear as a function of frequency, bridging the Brownian and non-Brownian regimes, may 
provide insights for the understanding of systems with more complicated interparticle inter- 
actions such as pastes, slurries, particle gels, jammed emulsions and metallic glasses under 
a wide range of conditions. 

We thank A. B. Schofield for the particles and acknowledge funding from EU FP7- 
Infrastructures 'ESMl' (CP&CSA-2010-262348) and Greek "Thales" project "COVISCO". 
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